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Abstract. 

We refine and extend a programme initiated by one of the current authors [Science 
276 (1997) 88; Phys. Rev. D56 (1997) 7578] advocating the use of the classical energy 
conditions of general relativity in a cosmological setting to place very general bounds 
on various cosmological parameters. We show how the energy conditions can be used 
to bound the Hubble parameter H(z), Omega parameter fl(z), density p(z), distance 
d{z), and lookback time T(z) as (relatively) simple functions of the redshift z, present- 
epoch Hubble parameter Hq, and present-epoch Omega parameter Qq. We compare 
these results with related observations in the literature, and confront the bounds with 
the recent supernova data. 

Dated: 10 December 2007; M^X-cd 2 February 2008 



Cosmodynamics: Energy conditions, Hubble bounds, density bounds, 



2 



1. Introduction 



Some 10 years ago one of the current authors initiated a programme of using the 
classical energy conditions of general relativity to place very general bounds on various 
cosmological parameters [H El |3]. In that early work, attention was mainly focussed 
on the energy density p(z) and lookback time T(z). Since then, the classical energy 
conditions have (on the one hand) seen continued use in studying issues such as 
the minimal requirements for cosmological bounces [U [5] and other "cosmological 
milestones" [6], and (on the other hand) have seen further applications to bounding 
cosmological distances d(z) [7J [8], and lookback time T(z) [9]. In the current article 
we shall try to draw these various threads together and establish several simple and 
rugged energy-condition-induced bounds on cosmological parameters. For some generic 
cosmological parameter, say represented by X(z), we shall seek bounds of the form 

X(z)^X hound = X f(Q ,z), (1) 

where Xq is the value of X(z) at the present epoch, the direction of the inequality may 
depend both on the bound being considered and the redshift region of interest, and 
/(f^O) z) is some dimensionless function to be determined. 

There is (at least) one important caveat: It should be kept clearly in mind that the 
classical energy conditions are not fundamental physics — in fact the classical energy 
conditions are known to be violated by quantum effects pUl EH EE2J [13], at least to 
some extent, and so the energy conditions should always be viewed provisionally — as 
a way of characterizing whether or not a certain situation is describable by "normal" 
physics [I21CE3]. 



2. Basic formulae 



In standard cosmology, one assumes the cosmological principle, that is, our universe 
is isotropic and homogeneous on large scales. This assumption leads one to consider 
cosmological spacetimes of the idealized FLRW form rTJJ [151 EE El [HI [19] : 

ds 2 = - dt 2 + a(t) 2 | j^-i + r 2 [d6 2 + sin 2 6 d0 2 ] j . (2) 

If we further assume that gravitational interactions at large scales are described by 
general relativity, we can use the Friedmann equations that relate the total density p 
and the total pressure p to a function of the scale factor a and its time derivatives. 
Indeed, in units where 8ttGn = 1, but explicitly retaining the speed of light c, we havdj] 

PW=3^ + ^), (3, 



P(t)= "2- -3--^, (4) 



a a 2 a 2 



X Note that we are now specifically assuming Friedmann dynamics for the universe, one is thus explicitly 
stepping outside the "cosmographic" or "cosmokinetic" framework of [201 12H 1221 1251 1231 1251 ?M \Tf[ [25] . 
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p(t) + 3p(t)= -6-. (5) 
a 

The classical energy conditions of general relativity, to the extent that one believes 
that they are a useful guide [121 allow one to deduce physical constraints on the 
behaviour of matter fields in strong gravitational fields or cosmological geometries. 
These conditions can most easily be stated in terms of the components of the stress 
energy tensor T^ v in an orthonormal frame. Ultimately, however, constraints on the 
stress-energy are converted, via the Einstein equations, to constraints on the spacetime 
geometry — in particular in a FLRW spacetime one is ultimately imposing conditions on 
the scale factor and its time derivatives. For a perfect fluid cosmology, and in terms of 
pressure and density, the so-called Null, Weak, Strong and Dominant energy conditions 
reduce to [29] : 

NEC: p + p > 0. 

In view of the Friedmann equations this then reduces to 

a a 2 he 2 , a a 2 he 2 

-- + ^ + ^>0; that is - < -= + 6 

WEC: This specializes to the NEC plus p > 0. 

This then reduces to the NEC plus the condition 

a 2 + kc 2 >0. (7) 

This condition is vacuous for k e {0, +1} and only for k = — 1 does it convey even 
a little information. 

SEC: This specializes to the NEC plus p + 3p > 0. 

This then reduces to the NEC plus the deceleration condition 

a , , 

- < 0. (8) 

DEC: p±p > 0. 

This reduces to the NEC plus the condition 

a (a 2 kc 2 \ 

Note particularly that the condition (jHJ) is independent of the space curvature k. Now, 
DEC implies WEC implies NEC, and SEC implies NEC, but otherwise the NEC, WEC, 
SEC, and DEC are mathematically independent assumptions. In particular, the SEC 
does not imply the WEC. Violating the NEC implies violating the DEC, SEC, and 
WEC as well [29]. 

Note that ideal relativistic fluids satisfy the DEC, and certainly all the known forms 
of normal matter encountered in our solar system satisfy the DEC. With sufficiently 
strong self-intereactions relativistic fluids can be made to violate the SEC (and DEC); 
but classical relativistic fluids always seem to satisfy the NEC. Most classical fields 
(apart from non-minimally coupled scalars) satisfy the NEC. Violating the NEC seems to 
require either quantum physics (which is unlikely to be a major contributor to the overall 
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cosmological evolution of the universe) or non-minimally coupled scalars (implying that 
one is effectively adopting some form of scalar-tensor gravity). 

Using this dynamical formulation of the energy conditions, Santos et al. [7j derive 
some bounds, for the special case k = 0, on the luminosity distance di of supernovae, 
and then contrast this with the legacy [301 El] and gold [32] datasets. In reference [Ej 
bounds on the distance modulus are presented for general values of k G {— 1,0, +1}, 
while in reference [9] they concentrate on the lookback time. Herein, we shall use a 
similar but distinct approach to obtain rugged bounds on the Hubble parameter, the 
Omega parameter, the density, the lookback time, and on the various distance scales 
defined previously in [261 EZ] for all values of k space curvature e {—1,0,+!}. 



3. Bounds on the Hubble parameter 



The energy conditions translate, in a FLRW setting, into the inequalities (jEJ), (JZJ), (jHJ), 
and (Q, from which we deduce bounds on the Hubble function H(z) in terms of the 
Hubble parameter Hq, the Omega parameter Qq, and the z-redshift. 

NEC: Using inequality we obtain: 



a da \a 



kc 2 



which can be integrated to yield 



That is 




da < 



kc 2 



da. 



Now using 



Hi 



a 



H{z) 2 < -kc 2 {a c 



l + z, 



and the relation 



1 + 



kc 2 



a H 



(10) 

(11) 

(12) 
(13) 

(14) 



after a few rearrangements we obtain the bound: 

H(z) > H NEC = H 0V /Q + [l-Q ](l + z) 2 . (15) 

In order to obtain this inequality we have assumed that z > 0, so that one is looking 
into the past. When looking into the future, z < 0, the inequality is reversed^ 

§ Note that looking back into the past z > 0, with z = oo corresponding to the big bang. In contrast, 
looking forward into the future z < 0, with z = — 1 corresponding to infinite expansion [201 127] . 
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Technical point: Existence of y^o + [1 — fi ] (1 + z) 2 in the K. domain. 
The expression 

fi + [l-^o](l + ^) 2 >0, (16) 

holds for all values of z, if and only if, Q ^ 1. In contrast, note that when Q > 1, 
there exists a z value for which the expression in the square root becomes negative 
or zero. This specific value of the z-redshift is given by: 



^ C = V^T" L (17) 

Note that at z = z NEC , we get H NEC (z NEC ) = 0. Also note that z NEC is positive 
as long as Qq is positive. Nothing unusual need happen to the universe itself at 
z NEC , it is only the bound that loses its predictive usefulness. In practice, given that 
current observational estimates are 

fl = 1.02 ±0.02 (PDG 2004 [33]), (18) 
n = 1.003+g;^ (PDG 2006 [3|), (19) 

let us point out that z NEC (Q = 1.04) = 4.1 and z NEC (Q = 1.01) = 9.0. So given 
current observational estimates of Qq, the fact that the highest-z supernovae seen 
to date have z < 2, and the fact that we are expected to run out of galaxies by 
the time we reach z < 7, the limitations associated with z NEC are unlikely to be 
significant in any realistic setting. 

Overall, the bound on the Hubble function (fl5j) is valid for Q < 1, V z G [0, +oo), 
and for Qq > 1 under the condition that z G [0, z NEC ]. 

WEC: From inequality ([7]), we can deduce that 

for k = — 1, a < V^k c. (20) 

To obtain a relation on the Hubble function, we divide equation (|20|) by a, and 
obtain the rather weak bound: 

H(z) > H WEC = H (1 + z) Vl - Q VO G(0,l). (21) 

We have assumed that z > in inequality (|2ip so that one is looking into the past. 
When looking into the future z < 0, the inequality is reversed. Note that this 
bound is only valid V fl G (0, 1) and Vz>0. 

Important remark: Note that as long as f2 > we have: 

H NEC > H WEC . (22) 

That is, the WEC really does not give us anything extra beyond the statement that 
Qo is positive. 

SEC: From inequality (|HJ), we deduce that 

1 1 

Va < a - < — . (23) 

a H a 
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Further, to obtain a relation on the Hubble function and Hubble parameter, we 
multiply equation (123]) by a, and we obtain the bound: 

H(z)>H SEC =H (l + z). (24) 

We have assumed that z > in inequality (I24I) so that one is looking into the past. 
When looking into the future z < 0, the inequality is reversed. Note that this 
bound can also be found in [3, [8] . 

DEC: To satisfy this energy condition, the NEC must hold as well as inequality (J9j). 
We use the same approach as for the NEC, rewriting (Q as: 

d(a 2 d) 



dt 

that is, 

da d 



2kc 2 a > 0, (25) 



v a 2 a) + 2kc 2 a > 0. (26) 
dt da 



Multiplying by a , this inequality leads to 



d_ 

da 



2 



1 / 2 -\2 . kc 

- (a aj H — — i 



> V a. (27) 



2 v ' 2 
Integrating, we can deduce the new inequality, 

V a < a (a 2 a) 2 + kc 2 a 4 < (a^ao) 2 + k c 2 a^. (28) 

Now, we multiply or divide appropriately by some combination of a and a to force 
the appearance of the Hubble function H(z) and the Hubble parameter Hq. We 
also use equations (Tl3l) and (JHJ) and substitute, leading to: 



H(z) < H DEC = Ho (1 + zWl + n [(l + z) 4 -l], (29) 

Again, we have assumed that z > in inequality fl29l) so that one is looking into 
the past. When looking into the future z < 0, the inequality is reversed. 
Thus the DEC is satisfied if and only if: 

H NEC < H(z) < H DEC , (30) 

where H NEC and H DEC are defined respectively in equations (ESP, and (GSP- 



Technical point: Existence of y/l + fio [(1 + z) A — 1] in the K domain: 
The expression 

1 + [(l + ^) 4 -l] >0, (31) 

holds, for all values of z > — 1 (and hence z > 0) if Qq G (0, 1). However, note 
that when ^ 1, there exists a 2 value for which the expression in the square root 
becomes negative or null. This specific value of the z-redshift is given by: 

o i\l/4 
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Note that z DEC is always negative so it is never a problem when looking back into 
the past. In fact z DEC (fl = 1.04) = —0.56 and z DEC (f2 = 1.01) = —0.68 are well 
into the future. 

Thus the bound on the Hubble function (129]) is valid for Q G (0, 1), V z > —1, and 
for Q ^ 1 under the condition that z G [z DEC , +oo]. If we are only interested in 
looking into the past, then the DEC bound holds for f2 > and z > 0. 



4. Bounds on the distance scales 



In order to obtain bounds on the various distance scales, it is enough to obtain a bound 
on Peebles' angular diameter distance [T7j and then use the different relations between 
the various distance scales presented in [26l 127] fjjl We choose to work primarily with 
Peebles' angular diameter distance because it minimizes the number of factors of (1 + z) 
occurring in the various formulae. Peebles' angular diameter distance is generally defined 
in its exact form as [D 



where 



d P (z) = a sin k { [ —A- dz } , (33) 

H a J H{z) 



sin(x), k = +1; 

sin fc (x) — \ x, k = 0; (34) 

sinh(x), k = —1. 

By changing variables and adopting definitions as in equations ( |T3l) and (I14j) . we 
can rewrite Peebles' angular diameter distance in an alternative exact general form, 
V z G [-1, +oo) and V fixed fi 8 

sinh [y/T=W f* ^ dz 



dp{z 
where we note 



c 



Hq \/l — O 



(35) 




k = +1; 

n { =1, k = 0; (36) 
k = -1. 

Observe that by continuity of the functions sinx/x and sinhx/x as x — >• 0, the function 
dp(^) is also continuous as Q — > 1 ± . For convenience, from equation (1331) . the angular 
diameter distance is given by 

, , , c sinh [VI - J] 

MZ) = JT y/l-tto ' (37) 

|| Peebles' angular diameter distance is equal to Weinberg's proper motion distance [15], and is also 
equal to D'Inverno's version of luminosity distance [35| . Details on how the various distance scales are 
inter-related can be found in [26l [27] . 

Another notation that is sometimes used is fi^ = 1 — f2o, so that k = — sign(f2fc). 
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where J is the integral defined by 

t f z Hq pda 

J= / — — dz = H a / — . (38) 
Jo H {z) J a a a 

The procedure now is as follows: The energy conditions provide bounds on H(z), which 

allow us to obtain a bound on the integral J. Then provided the function sin^ is 

monotonic on the interval z G [0, +oo], (or at least some sub-interval z G [0, z max ] ), we 

can derive a bound on the angular diameter distance on this same domain. 



NEC: The null energy condition gives a bound on H in equation (1151) leading to the 
inequality, 

J= / tttM* < J NEC = / ; (39) 
J H(z) - J v/fio + [1 - fi ] (1 + zf 

We integrate, and substitute the resulting bound back into the angular diameter 

distance. In the general case, we obtain the bound 



d P {z) < d PNEc = [l + z- ^n + (i-n )(i + zy\ ; (40) 

Vfi < 1, Vz G [0,+oo]; 
VOo>l,Vze [0,2 NBC ). 

Note that as f2o — > 1 (A; = 0), we have 

dp(z)<d LNEc = ^- n = l; VzG[0,+oo]. (41) 

so we find the same particular result as in [7], that is, 

d L {z)<d L ^ c = CZ[l H + Z) no = l; V*e[0,+oo]. (42) 

Equation (1401) is the more general case, valid for all values of k G {— 1,0, +1}. 
(Equation (15) of reference j8] can be viewed as an intermediate step in deriving 
equation (JIT?]) above.) 

Note that as Qq — * 1, equation (140]) can be developed in a Taylor series as 

dp^iz) = Jf + ^- Q ^-l) + O ([flo - l] 2 ) . (43) 
If instead one performs a low-redshift expansion, then for general n 

cz f 1 (n — 1) 2 



^ NBC (*) = ^ ( 1 + + ^) ) • (44) 

WEC: The weak energy condition gives a (weak) bound on H(z) in equation (l2~Tj) only 
for n G (0, 1) leading to the inequality, 

We integrate, and substitute the resulting bound back into the angular diameter 
distance. We obtain 

d P {z) < dp = p Z f ^ + Z } \ (46) 

WEC H VT^TTo (1 + *) 



vn g (o,i), v z e [o,+oo]. 
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Important remark: Note that 

dp <d P . (47) 

NEC — WEC v ' 

Thus the bound dp is not very useful. 

WEC J 

SEC: This energy condition gives a bound on H(z) in (J24l) . and therefore 

J<J SEC = / _ = ln(l + z). (48) 

JO 1 + z 

In the general case, we obtain the bound on the angular diameter distance (cf. 
equation (17) of [5]): 

dp(z) < dp w = ±. si "Mv / ™^d + 4 ; (49) 

VO < l,Vz G [0,+oo]; 
Vfi > 1, Vz G [0 j ^max I • 

In particular, 

• For k = -1, that is V fi < 1, and V z E [0, +oo], 

dp ^ w - w uf^m ' ( } 

• For k — 0, that is when Qq = 1, and V z G [0, +oo] (cf. the equivalent result 
in H), 

dp ^ {z) = ir ln ^ 1+;z ) ; ( 51 ) 



• For fc = +1, that is V f2o > 1, and V z G [0, -2 max ], 

c sin [y^T In (1 + zj\ 

dp ^ {z) -w Tn^r ' (52) 

whereE] 

W = exp(^ = =)-1. (53) 

In fact, since z max (fl = 1.04) = 2575 and z mabK (Q = 1.01) = 6.6 x 10 6 , we see 
that this constraint is never a significant limitation on the bounds. 

Note that as Qo — > 1 + , z max — > +oo, and equation ( 1491) can be developed in a Taylor 
series as 

d PsEC (z) = ^- ln(l + *)-^L [ln(l + ,)] 3 (fi -l) + O([fi -l] 2 ). (54) 



If instead one performs a low-redshift expansion, then for general Qq 

dp ^) = wA l - Z 2 + °^ 2) }- (55) 



+ If J < J SEC then for the sine function we have sin(\/f2o — 1 ^) ^ sin(\/f2o — 1 ^sec)) provided that 
< \/VL — 1 J SEC < 7r/2, leading to the condition that 2 < z max = exp ( 2% /q -i ) — ^' 
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DEC: Remember that to satisfy the DEC, the Hubble function needs to satisfy both the 
NEC, inequality (fl5j) . and the second inequality (j29j) . As a consequence, in order 
for the DEC to hold, Peebles' angular diameter distance must satisfy inequality 
(140|) . and a second inequality to be derived below. From equation (1291) . we obtain 



/ ^ < J DEC ^ / = ■ (56) 



/„ H(z) - DEC J„ (l + z)^l + n„ \(l + z)'-l] 

This integration is a bit more tricky than the previous integrations for the NEC 
and SEC0 We obtain: 

^ * ,J d-".*^ (! + «)' y (57) 
2V1 -ilo [i-n + VT^TT Jl + fi [(1 + zf - 1] 



In the general case, we obtain the lower bound on the angular diameter distance, 



dpiz) > dp (z) = 



DEC 



1 + O [(1 + zf - 1] - (1 + fi [(1 + zf - 1]) 



#o(l + ^n 2fi (1 - 



VO < 1, M z G [0,+oo]; 
VJlo>l,Vze (z DEC ,+oo]. 

(Equations (19) and (20) of [8] can be viewed as intermediate stages in deriving 
this result.) Note that contrast to the situation for the SEC, there is no constraint 
on a maximum value for z coming from the requirement that the sine function be 
monotonic. 

The lower bound of the DEC in equation (|58j) can also be represented in a Taylor 
series as fio — * 1, 

d (z\- C Z{2 + Z) l ^ 2 (2 + ^) 2 (3^ + 6^ + 4) n , orro m 

(59) 

If instead one performs a low-redshift expansion, then for general Qq 

^,„W = f{l- e ^ + ^)}. (60) 



Energy conditions and Supernovae data: We can plot the angular diameter distance 
bounds (NEC, WEC, SEC, and DEC) and compare them with the angular diameter 
distance raw data from the supernova datasets. We have used data from the supernova 
legacy survey (Iegacy05) [30],[3l] and the Riess et. al. "gold" dataset of 2006 (gold06) [32J. 

• Figure [1] compares the upper bounds (NEC and SEC), and the lower bound (DEC), 
with the Iegacy05 dataset. In contrast figure [2] uses the gold06 dataset. 

* Indeed, common symbolic manipulation systems such as Maple or Mathematica require significant 
human intervention before they will even condescend to verify these results. 




Figure 1. This figure shows Peebles' angular diameter distance dp(z) as a function 
of the z-redshift from the nearby and legacy survey Iegacy05 dataset [30], [31] . Data 
under the "red solid lines" satisfy the NEC/WEC, data under the "black dashdot lines" 
satisfy the SEC, data under the "red solid lines" and above the "magenta dashed lines" 
satisfy the DEC. The 5 lines for each energy conditions correspond to varying values 
of the parameter fi = {0.96, 0.98, 1.00, 1.02, 1.04}. The value of the Hubble constant 
is taken to be Hq = 70 km/s/Mpc. 



• To satisfy the NEC, the data must be under the "red solid" bound, and we can see 
that most of the data seem to satisfy this condition. 

• For the SEC to hold, the data must be under the "black dashdot" bound. Visually 
it seems "obvious" that the data significantly violate the SEC. 

• Finally, the DEC is satisfied if both: (1) the NEC is satisfied, and (2) if the data 
is above the "magenta dashed" lower bound. This latter condition is well satisfied 
for the bulk of the data, therefore satisfying the DEC is dependent on the NEC 
holding. 

• As is traditional for estimates of cosmological distance, we plot only one-sigma 
statistical uncertainties, without any allowance for systematic uncertainties. Any 
realistic attempt at more careful treatment of the systematics, and/or going to 
3-sigma error bars, makes the plots much more disquieting [26] [27] . 

• There seem to be noticeable visual differences when looking at figures [1] or [2] which 




Figure 2. This figure shows Peebles' angular diameter distance dp(z) as a function 
of the z-rcdshift from the gold06 dataset [35]. Data under the "red solid lines" satisfy 
the NEC/WEC, data under the "black dashdot lines" satisfy the SEC, data under 
the "red solid lines" and above the "magenta dashed lines" satisfy the DEC. The 
5 lines for each energy conditions correspond to varying values of the parameter 
fl Q = {0.96,0.98,1.00,1.02,1.04}. The value of the Hubble constant is taken to be 
H = 70 km/s/Mpc. 



make it tricky to conclude whether the classical energy conditions are satisfied or 
not by just looking at the supernova data in isolation. For example, there are a few 
supernovae data in the redshift range 0.8 < z < 1 that appear to violate the NEC 
in an obvious manner for the Iegacy05 dataset in figure ([1]) . However, the violation 
does not appear to be as dramatic when looking at the same range of data in the 
gold06 dataset in figure ([2]). Another example is that the NEC naively seems to be 
violated for data in the redshift range 0.4 < z < 0.6 in the gold06 dataset. On the 
contrary, even if there are less data in the Iegacy05 dataset, one cannot draw the 
same conclusion. 

• Unlike references [71 [8], [9] , we believe we cannot draw any conclusions by using the 
low-redshift linear part of the distance scale curve, as the data has been scaled to 
enforce a particular value of H . 

• Note in particular that for low redshift the luminosity distance is bounded, both 
above and below, by constraints of the form cz/H + 0(z 2 ). Thus for data with any 
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statistical uncertainties whatsoever, at low enough z, one would expect roughly half 
the supernovae to "violate" one or more of these bounds. 

• It is important to realise that the "slope" of the bounds at z = depends on the 
estimate of Hq one adopts (and the precise value of c) , and can further be affected 
by the value of the magnitude "offset" reported for the data. 



5. Bounds on the lookback time 



The lookback time is defined as [TJ [2j [3] 



T(z) 



That is 



T(z) 



f dt , 
dt= — da 
J da 

1_ d[a /(l + z)] 

H a /(l + z) 

1 



a da 
a a 

1_ dz/(l + z) 2 
H 1/(1 + z) ' 



dz. 



(61) 



(62) 



(1 + z) H{z) 

In order to obtain bounds on the lookback time for the different energy conditions NEC, 
WEC, SEC, and DEC, we again use the bounds on the Hubble parameter H(z). 

NEC: The null energy condition gives a bound on H(z) in equation (|15[) . leading to 
the inequality 



T{z) 



1 + z) H(z) dZ ~ T ^ {Z) 



# NEC w 



dz 



(63) 



We integrate, and substitute the resulting bound into the lookback time (equivalent 
formulae can be found in reference [2], see also equation (14) of [9]) 



T(z) < T NEC (z) 



111 



;i + z) (i + vn~ ) 



(64) 



this bound is valid V Q < 1, V z > and V Q > 1, V z e [0, z NEC \ 
Alternatively, this bound can be rewritten as 



T(z) < T N 



1 



H \/Q, 



In 



o 



n + (i-n )(i + z f- 

(1 - y/H^ (1 + Z) 



Using the standard result that sinh 1 x = ln(x + \/ x 2 + 1), we can for k 
is, Q < 1) also re-cast this as 



(65) 
-1 (that 



T{z) < T N 



sinh 1 



— sinh 1 



(66) 
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Equivalent formulae can be found in reference [2]. Similarly for k — +1 (that is, 
f2o > 1) we can use the fact that cosh" 1 x = ln(x + \J x 1 — 1) to obtain 



T(*) < T NEC (.) = <[ cosh" 1 I - cosh- 1 / 



'(67) 

(Note that this formula is valid only for z < z NEC , since otherwise the argument 
of the cosh -1 is less than unity.) Finally, the upper bound derived from the 
NEC in equation (1641) . or equivalently any of equations (!65|) -(!66 |) - (l67|) . can also 
be represented in a Taylor series as Qq — > 1: 

W*) = + ^ + 2,-2ln(l + ,) _ flo _ 

-no 4ii 
In particular for f2 = 1 we recover the results of [2] and [9]. 

WEC: This energy conditions gives a bound on #(2) in f[2"Tj) for f2 £ (0, 1) only 
Thereby it can be deduced that 



T{z) = L (i + 7) W) iz - T ^ (z) s I (i + 7) O (69) 

providing the (weak) bound 



T(*) < T WEC ( 2 ) = — — -====————. (70) 
#0 Vl - fio (! + 2 ) 

Again, this provided no additional useful information beyond the NEC-derived 
bound. 

SEC: This energy condition gives a bound on H(z) in (124"|) . Thereby it can be deduced 
(as in the articles P El El EJ) that, 

f* 1 r 1 

Jo {l + z)H{z) Jo ^ + z)H SEC {z) 

that is 



r -° w = TmTTI)- (72) 

The above result equation (I72j) is completely independent of the value of the 
parameter Qq- Note that this result was first introduced by Visser in [U [2j [3]. 

DEC: Remember that to satisfy the DEC, the Hubble function needs to satisfy the 
NEC, inequality (fT5j) . and inequality (I29l) . As a consequence, in order for the DEC 
to hold, the lookback time must satisfy inequality (I64j) . and a second inequality 
that we shall derive below. From equation (1291) . we obtain 

T(2) = I (i + ') m dz ~ Tdeo(2) s I (rnyiww dz> (73) 

that is 

1 /" z 1 
T DEC (z) = — / d*, (74) 

DECV H J (1 + ^)2 ^1 + f2 ((1 + 2?) 4 - 1) 



In 



(76) 
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The integration of this bound is considerably harder than for the other energy 
conditions. Again, common symbolic manipulation systems such as Maple or 
Mathematica require significant human intervention before producing anything 
useful. Let us first write 

1 [* 1 
T„„„m = -== I . dz, (75) 

DEcV ' HoV^oJo (1 + zY v /l-(l-^o 1 )(l + ^)- 4 
and then, (following the procedure of [2]), apply the binomial theorem 

[i - (i - iv) (i + z)- 4 ] - 1/2 = f; (-y 2 ) (-it (i - n^r (i + z) 

n=0 ^ ' 

Now the binomial series will converge provided 

\(l-Q^)(l + zr*\<l : (77) 

and in view of the region we are integrating over, this means that the series for the 
lookback time will converge provided 

1 1 — x | < 1, that is fi e(l/2,oo). (78) 

Subject to this condition we can integrate, and obtain the convergent series 

As a practical matter, for many purposes this series representation may be enough, 
but we can tidy things up somewhat by first defining 



_ n V n J 4n + 3 



(80) 



n=0 

in which case 



(l + z)< 



Finally we can recognize that S(x) is a particular example of hypergeometric serie 
and so write 

Therefore 



T DEC (z)= (84) 
(t The classical hypergeometric series is given by 

tp f 7 \ ( a )n(^)n X , , 

2J F 1 (a,6;c;x) = ^— — , (82) 

n=0 lCj " 

where (a) n = a(a + l)(a + 2) (a + n — 1) is the rising factorial, or Pochhammer symbol. The series is in 
general a convergent power series for values of x such that \x\ < 1. 
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Again this agrees the results reported in [2J, and more recently (for k — 0, flo — 1) 
in [9]. Of course writing the result in terms of hypergeometric functions does 
not necessarily give one much additional physical insight — for physical insight 
the series S(x) is sufficient, and the realization that one is in fact dealing with a 
hypergeometric function is likely to be useful only if for some reason one wishes 
to numerically programme the bound into a computer. (The result can also be 
cast in terms of elliptic integrals [2J, but this does not appear to be particularly 
illuminating.) 

6. Bounds on the Omega parameter Q(z) 
We have the following identity 

U 2 U r 2 2 tt2 tt2 

That is 

fl(z) = 1 + (flo - 1) {l + zf j^L. (86) 

Therefore, a bound on H(z) automatically implies a bound on Q(z). 

NEC: The null energy condition gives a bound on H(z), as in equation (1151) . leading 
to 

" NEC = fio + (l-^)(l + ^) 2 ' (87) 
and the inequalities 

ifft o <l,V*>0, thenQ(z) > Q NEC , (88) 

iffi = l,V2>0, thenfi(^) = Q NEC = 1; (89) 

iffi >l,V2>0, thenfi(z) < fi NEC . (90) 

Note that as £Iq — > 1, equation (1ST)) can be developed in a Taylor series as 

fi NEC = 1 + (1 + zf (flo - 1) + O ([flo - 1] 2 ) • (91) 
WEC: The weak energy condition gives a bound on H(z), as in equation (I2TT) . but for 

no G (0, 1) only, leading to the trivial result n wEC = 0, and the trivial inequality 

iffi <i,v*>o, n > n WEC = o. (92) 

This bound is not useful, except as a consistency check. 
SEC: The strong energy condition gives a bound on H(z), as in equation (124)1 . leading 
to 

n sEC = n , (93) 
and the inequalities 

if n < 1, V z > 0, then Q(z) > Q SEC = Q ] (94) 
ifn = l,V2>0, then n(^) = n sEC = n = 1; (95) 
if n > 1, V z > 0, then Q(z) < Q SEC = Q - (96) 
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DEC: The dominant energy condition gives a bound on H(z), as in equation 
leading to 



1 + O0 [(l + zf-iy 

and the inequalities 

iffi <l,Vz>0, then Q(z) < Q DEC ; 

if Q = 1, V z > 0, then = fi DEC = 1; 

if fi > 1,V z G [0,z DEC ], then Q(z) > fi DEC . 
Note that as fio — > 1, equation ([57]) can be developed in a Taylor series as 

1 , (^0-1) 



O ([Oo - l] 2 ) . 



(97) 

(98) 
(99) 
(100) 

(101) 



(1 + -) 4 

These bounds on Q(z) are potentially of interest with regard to cosmological 
nucleosynthesis, which is effectively sensitive to ^^nucleosynthesis)- More generally, any 
bound on the number of relativistic particle species at any particular epoch can be 
converted, with a little work and some technical assumptions, into a bound on the 
Omega parameter at that epoch. 



7. Bounds on the density p(z) 

We have the following identity 



P 



3QH 2 



lHl + z 



2 Hi 

H 2 



That is 

p{z) = 3H(z) 2 + 3 (ft - 1) (1 + ^) 2 # 2 , 

showing that a bound on H(z) automatically implies a bound on p(z). 
Alternatively, we can also write the following identity 

a 2 J ~ d # 2 + a 2 a 2 " fi JJg + Po ^ fi 

That is 

J_ i/(^ 2 A 1 



p = 3 ( if 2 + 



pO) = Po 



if 2 



1 



Again, a bound on ff(£) automatically implies a bound on p(z). 

NEC: The null energy condition gives a bound on H(z), as in equation ([TBI) 



to 



Pnec - 3 ^o Hq 



Po, 



and the inequality, 

V fio, V z > 0, p(z) > p NEC = po- 
This inequality was also derived by more direct means in [TJ [2] . 



(102) 
(103) 

(104) 

(105) 

, leading 
(106) 

(107) 
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WEC: The weak energy condition gives a bound on H(z), as in equation (12 ip . leading 
to 

Pwec = 0, (108) 
and the inequality, 

Vfi„<l,Vz>0, p > p WEC = 0. (109) 

Of course, since by assuming the WEC we have already assumed that p > 0, this 
bound is not very useful (and is at best a consistency check on the formalism). 
SEC: The strong energy condition gives a bound on H(z), as in equation f)24p . leading 
to 

p SEC = 3Q Hi (1 + zf = po (1 + zf , (110) 
and the inequality, 

Vfio.VoO, p(z) > p SEC = po (1 + zf. (HI) 
This inequality was also derived by more direct means in [U [2] . 
DEC: The dominant energy condition gives a bound on H(z) in equation (|29|) . leading 
to 

p DEC = 3Q Hi (1 + zf = po (1 + zf , (112) 
and the inequality 

Vn ,Vo0,p< p DEC = po (1 + zf . (113) 

This inequality was also derived by more direct means in (TJ [2] . 

Note that bounds on the density and the Hubble function are intimately related. A 
bound on one will automatically provide a bound on the other, and comments made 
regarding the Hubble bounds can be carried over to this situation as well. 

8. Bounds on the pressure p(z) 

For the pressure p(z) things are a little different; we have the following identity involving 
the second time derivative of the scale factor: 

a 2 kc 2 a 

p=--- — -2-. 114 



or a a 



But 



a 2 | 2 d _ 1 d(a 2 a) _ 1 d(a 2 a) _ 1 d[H 2 a 3 ] 
a 2 a da 2 dt a 2 da a 2 da ' 



implying 

Here the point is that one would need a bound on the derivative of H(z) in order to get 
a direct bound on the pressure p(z). This does not appear to lead to anything useful. 

However, if one has a bound on H(z) and hence p(z), one can indirectly get a 
constraint on p(z) via the classical energy conditions. Again, this does not appear to 
lead to anything useful. 
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9. Conclusions 



In this article we have extended the discussion of [U [2, [3], and more recently of [3, El E], 
to develop a number of rugged and general energy-condition-induced bounds on various 
cosmological parameters, bounds which have all taken the form 

X(z)^X hound = X f(n ,z) J (117) 

where X(z) is some cosmological parameter, X is its present-day value, and f(£lo, z) is 
a dimensionless function depending on the particular bound under consideration. While 
most of the bounds we have considered can be derived by "elementary" means, in at least 
two cases the integrals (and other algebraic manipulations) were sufficiently tricky that 
common symbolic manipulation programmes required considerable human intervention 
in order to obtain useful results. 

Overall, we have now managed to place bounds of the above form on the 
Hubble parameter, the Omega parameter, the density, the lookback time, and various 
cosmological distances - - where for simplicity of presentation we have focussed on 
Peebles' angular diameter distance as being representative. We have briefly sketched 
how to use these bounds to confront the supernova data, but have not yet performed 
any detailed analysis of this point. 
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